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Final  report  on  ONR  contract  #N00014-84-C-0261 


Summary  of  work  completed  under  the  contract 

During  the  course  of  this  contract,  research  has  been  conducted  in  three  distinct 
areas,  all  involving  the  Josephson  effect: 

A.  Development  of  Josephson  junction  arrays  as  tunable  sources  of  microwave 
through  submillimeter  radiation,  *, 

B.  Conformation  of  general  relativistic  effects,  equivalent  to  the  red  shift,  for 

/* 

charged  particles,  / 

r'”  “*  ... 

G.  Studies  of  the  quantum  mechanics  of  macroscopic  variables,  in  particular 
of  the  flux  linking  a  SQUID.  - _ — 

The  fundamental  studies  of  Josephson  oscillators  A.),  which  had  been  supported 
under  previous  ONR  contracts,  were  concluded  during  the  first  part  of  the  contract. 
The  present  status  of  this  work  is  discussed  in  detail  in  the  reprint  attached  as  reference 
A.  The  success  of  this  work  has  lead  to  very  active  ongoing  research  programs  in  my 
laboratory,  supported  by  several  sources  to  develop  technical  applications  of  these 
oscillators.  In  particular,  techniques  for  coupling  junctions  distributed  over  many 
wavelengths  were  developed  and  applied  to  the  design  of  a  broadly  tunable  submilli¬ 
meter  source  generating  7/iW  of  power  at  400GHz.  In  addition,  recent  results  on 
frequency  agile  microwave  sources  demonstrate  modulation  of  a  16GHz  carrier  over 
3GHz  at  modulation  frequencies  in  excess  of  500MHz. 

Our  focus  during  the  middle  years  of  the  contract  was  on  B.),  the  study  of  effects 
predicted  by  general  relativity  on  charged  particles  moving  through  a  gravitational  field. 
These  experiments,  which  are  described  in  the  reprint  in  appendix  B,  produced  the  first 
demonstration  that  the  electrochemical  potential  difference  between  two 
superconducting  wires  varies  with  gravitational  potential  (by  about  1  part  in  1016  per 
meter  in  the  earth’s  field).  This  effect  is  essentially  the  equivalent  for  charged  particles 
of  the  well  known  gravitational  red  shift  for  photons.  Use  was  made  of  essentially  ideal 
Josephson  effect  “batteries”  phase-locked  to  a  microwave  source  to  establish  precisely 
related  EMFs  at  different  heights.  As  part  of  this  work  techniques  were  developed  to 
measure  potential  differences  of  10-22V. 

During  the  final  phase  of  the  contract,  we  began  work  to  test  predictions  that 
macroscopic  variables,  such  as  the  flux  through  a  SQUID  loop,  display  quantum 


mechanical  properties  such  as  tunneling  and  coherent  superposition  of  the  flux  wave 
functions  belonging  to  different  states.  This  work  has  been  continued  under  a  subse¬ 
quent  ONR  grant  and  her.  produced  the  first  evidence  that  quantum  theory  for 
transitions  between  the  fluxoid  states  in  a  degenerate  two  level  SQUID  is  correct.  The 
results  under  this  contract  are  presented  in  the  reprint  in  appendix  C.  Here  a  novel 
technique  is  demonstrated  for  using  a  two  junction  composite  SQUID  to  make  what  is 
essentially  a  “man  made  atom”  with  energy  levels  and  lifetime  which  can  be  controlled 
in  situ  by  the  experimenter. 
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This  gives  a  mixing  current  0  v-ve  0  v/2  so 

=  a^e  COS  <50,  (4.6)  P'S-  ®;  (a)  Junciioi)  /-F  curve  in  neighborhood  of  a  radiation-induced  step,  (b)  Phase-lockii 

Where  <50  s  ©  -  <yer.  ®'°n8-  S'eP’  Note  greatest  stability  is  for  <5©  =  n/2.  The  regi. 
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1.  JOSEPHSON  ARRAYS  AS  HIGH  FREQUENCY  SOURCES 
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Test  for  Relativistic  Gravitational  Effects  on  Charged  Particles 

A.  K.  Jain,  J.  E.  Lukens,  and  J.-S.  Tsai(a) 

Department  of  Physics,  State  University  of  New  York,  Stony  Brook,  New  York  1 1794 

(Received  14  November  1986) 


Experimental  results  are  presented  which  provide  the  first  measurement  of  the  effects  of  a  gravitation¬ 
al  field  on  charged  particles,  equivalent  to  the  red  shift  for  photons.  Two  Josephson-effect  batteries 
(K—  300  pV)  having  a  vertical  separation  .of  7.2  cm  are  connected  in  opposition  by  superconducting 
wires.  A  voltage  difference  of  2.35xl0-21  V  is  maintained  between  these  batteries  by  means  of  the 
gravitational  red  shift.  The  emf  around  this  loop  is,  however,  measured  to  be  less  than  1  x  10  ~22  V,  con¬ 
sistent  with  the  predicted  invariance  of  the  gravito-electrochemical  potential  along  the  wires. 

PACS  numbers:  04.80.+Z,  74.50.+r 


The  strong  equivalence  principle,  which  states  that  lo¬ 
cally  the  effects  of  a  uniform  gravitational  field  are  indis¬ 
tinguishable  from  those  due  to  an  accelerating  reference 
frame,  is  a  fundamental  postulate  of  general  relativity. 
A  consequence  of  this  equivalence,  along  with  special 
relativity,  is  that  the  rates  of  identical  clocks  running  in 
different  gravitational  potentials  will  differ.  This  varia¬ 
tion  in  clock  rates  implies  the  gravitational  red  shift  of 
photons,  which  was  first  observed  by  Pound  and  Rebka 1 
and  has,  during  the  past  two  decades,  been  verified  by  a 
number  of  workers2  to  a  present  accuracy  of  better  than 
1  part  in  104.  There  have,  however,  previously  been  no 
experimental  tests  of  the  corresponding  effects  for 
charged  particles. 

One  such  effect  occurs  for  the  conduction  electrons  in 
a  metal,  such  as  a  superconductor,  in  a  gravitational  field 
where,  in  equilibrium,  the  electrochemical  potential  ft 
will  in  general  not  be  constant.  Rather,  a  related  quanti¬ 
ty  {i=fi(  1+1)  called  the  gravito-electrochemical  poten¬ 
tial  will  be  constant  along  the  wire.3'6  Here  X=zg/c 2 
and  depends  on  the  gravitational  potential  through  the 
height  z;  g  is  the  gravitational  acceleration  at  the  earth’s 
surface  and  c  the  speed  of  light  in  a  vacuum,  giving 
X/z  =  1.09xl0“16/m.  For  a  resistanceless  metal  such  as 
a  superconductor,  p  should  also  be  independent  of  posi¬ 
tion  even  with  a  dc  current  flowing.6  Thus  a  circuit  with 
two  identical  batteries  (as  measured  by  local  observers) 
of  emf  V,  separated  in  height  by  z  and  connected  in  op¬ 
position  by  superconducting  wires,  would  have  a  net  emf 
AV=XV  around  the  zero-resistance  loop  leading  to  a 
loop  current  increasing  linearly  in  time.4 

Within  a  single  wire  a  purely  Newtonian  gravitational 
field  can  produce  a  number  of  effects  which  shift  the 
electrochemical  potential.  For  instance,  the  gravitational 
force  on  the  electrons  must  be  balanced  by  an  electro¬ 
static  force  giving  rise  to  an  electric  field  (the  Schiff- 
Barnhiil  field7).  This  latter  effect  is  in  fact  contained  in 
the  potential  ft  if  the  electron  rest  mass  is  included,  since 
Xtnc 2  =mgz.  An  electric  field  can  also  be  produced  be¬ 
cause  of  the  shift  in  the  Fermi  level  induced  by  lattice 
distortion  due  to  the  gravitational  force  on  the  ions.s 


These  gravitational  effects  can  be  10  to  15  orders  of 
magnitude  greater9  than  the  relativistic  shift  in  potential 
difference  discussed  above.  However,  they  do  not  con¬ 
tribute  to  the  net  emf  around  a  loop. 

In  this  Letter  we  report  the  results10  of  a  null  experi¬ 
ment  to  compare  the  fractional  change  with  height  of  the 
electrochemical  potential  difference  between  two  super¬ 
conducting  wires  to  the  fractional  change  in  photon  fre¬ 
quency  over  the  same  height,  which  is  given  by  Av/v=X. 
The  schematic  of  the  experiment  is  shown  in  Fig.  1.  The 
“batteries”  used  to  provide  the  highly  precise  reference 
emFs  are  Josephson  junctions  separated  vertically  by  a 
height  z~1.2  cm  and  phase  locked  to  a  common  exter¬ 
nal  microwave  source.  The  relation  between  the  fre¬ 
quency  v  of  the  external  radiation  and  the  dc  electro¬ 
chemical  potential  difference  across  the  junction,  K=A/r, 
is  given  by  the  Josephson  equation  as 

K=nv<t>o.  (l) 

where  the  integer  n  is  the  order  of  the  radiation-induced 
step  in  the  junction  I-V  curve.  Since  the  radiation 
reaching  the  upper  junction  is  red  shifted  with  respect  to 


FIG.  1.  Schematic  of  the  measurement  circuit.  The  crosses 
denote  the  Josephson  batteries.  The  darker  lines  indicate  su¬ 
perconducting  portions  of  the  circuit. 


©  1987  The  American  Physical  Society 


1165 


Volume  58,  Number  12 


PHYSICAL  REVIEW  LETTERS 


23  March  1987 


that  reaching  the  lower  junction,  the  potential  differences 
across  the  upper  and  lower  junctions,  Vu  and  Vt  respec¬ 
tively,  will  differ  slightly  so  that 

K„  =  F/(1-U  (2) 

In  the  absence  of  additional  relativistic  effects  on  the 
Cooper  pairs,  one  would  thus  expect  the  net  emf  AV  in 
the  loop  containing  the  junctions  to  be  XV— 2x  10-2i  V, 
since  V— 300  pV  and  X— 7x10“ 18  for  this  experiment. 
Such  a  nonzero  value  for  AV  would  lead,  as  discussed 
below,  to  a  time-varying  flux  through  the  loop.  This 
could  be  detected  by  use  of  a  superconducting  quantum 
interference  device  (SQUID)  magnetometer,  as  shown  in 
Fig.  1.  The  relativistic  prediction  that  ff=const  in  the 
superconducting  leads,  however,  implies  that  the  poten¬ 
tial  difference  V  between  the  leads  varies  with  height  so 
that  K(z)  =  K(0)(l  —  Xq).  The  symbol  Xq,  which  in 
theory  equals  X,  has  been  introduced  to  allow  for  the  pos¬ 
sibility  of  corrections  to  the  predicted  relativistic  effects 
on  charged  particles.  The  loop  emf  is  thus  AV—V ) 
x (X  —  Xq)  and  is  predicted  to  be  zero  on  the  basis  of  the 
strong  equivalence  principle. 

The  basic  technique  for  the  detection  of  extremely 
small  voltage  differences  between  two  Josephson  junc¬ 
tions  by  monitoring  of  the  flux  change  in  a  superconduct¬ 
ing  loop  was  developed11  shortly  after  the  discovery  of 
the  Josephson  effect.  The  apparatus  used  for  these  mea¬ 
surements  is  a  refined  version  of  that  previously  used  by 
us  to  demonstrate  the  universality  of  Eq.  (1)  for  different 
types  of  Josephson  junctions12  to  2  parts  in  1016.  The 
Josephson  junctions  used  here  are  lead-alloy  tunnel  junc¬ 
tions  with  critical  currents  of  7C— 900  /iA.  The  junctions 
are  shunted  with  low-inductance  0.5-ft  resistors  to  en¬ 
sure  nonhysteretic  operation  and  have  I-V  characteris¬ 
tics  that  are  extremely  close  to  those  predicted  for  ideal 
resistively  shunted  junctions  (RSJ  model).  The  junc¬ 
tions  and  superconducting  leads  connecting  them  are 
placed  over  a  niobium  ground  plane  giving  a  total  loop 
inductance  of  £  =  1  nil,  associated  primarily  with  the 
coupling  inductor  to  the  magnetometer.  The  microwave 
radiation  is  coupled  to  the  junctions  by  our  placing  each 
junction  at  the  end  of  a  capacitively  grounded  microstrip 
which  is  in  turn  coupled  to  a  coaxial  cable  through  a 
small  section  of  coplanar  waveguide.  In  order  to  mini¬ 
mize  any  extraneous  shifts  in  the  relative  phases  of  the 
radiation  reaching  the  two  junctions,  these  two  coaxial 
cables  are  connected,  just  outside  the  sample  cell  and 
halfway  between  the  junctions,  to  a  common  cable  lead¬ 
ing  out  of  the  cryostat. 

In  general  Josephson  junctions  will  phase  lock  to  ap¬ 
plied  radiation  over  a  range  of  bias  current  A/,  “./(sin(<5, 
—  zr/2),  where  the  maximum  current  variation— or  lock¬ 
ing  range— is  dependent  on  microwave  power.  Here, 
8i=9j—n(Oit  is  the  difference  between  the 

linear  component  6,  of  the  phase  of  the  Josephson  oscil¬ 
lations  in  the  /th  junction  and  the  phase  of  the  applied 


radiation  at  the  junction. 13  At  a  constant  gravitational 
potential  the  time-averaged  fiuxoid  quantization  condi¬ 
tion  for  the  loop  of  inductance  L  as  shown  in  Fig.  1  is 

then14 


(Oo/2;r)(02-0, )+£//+<!><,  =w<l>o,  (3) 


where  d>a  is  the  flux  applied  to  the  loop  and  7/ 
=  (A7 2  —  A7 1  )/2.  For  simplicity,  identical  junctions  and 
locking  ranges  for  1  and  2  have  been  assumed. 

The  emf  which  drives  the  changes  in  the  loop  current 
7/  is,  from  Eq.  (3), 


$0  d  (s  c  \  .  A  , 

~z — —  (52“<5|)  +  —  nAco+ 
2 n  at  2k 


d<bg 
dt  ’ 


where  Aco=co2~ co\.  For  the  SQUID's  used  in  these 
measurements  P'i=2kLJ /Oo^  1,  thus  the  terms  in  8  can 
be  neglected,  giving 


.dlt  _  d<»a 
dt  dt 


, _®D  . 

+  — -  n  Aco. 
2k 


The  response  of  this  phase-locked  system  to  a  small  fre¬ 
quency  difference  Aco,  as  described  by  Eq.  (3),  is  the 
same  as  the  response  to  an  applied  flux  varying  linearly 
with  time  for  the  system  biased  below  the  junction  criti¬ 
cal  current  in  the  absence  of  radiation.  For  example,  if 
d>a  is  changed  (or  Aw^O),  a  loop  current  will  be  in¬ 
duced.  As  long  as  7/  <  J,  m  will  remain  constant  and  7/ 
will  increase  (nearly)  linearly  with  <J>a  (or  time)  until 
7 )—J,  then  m  will  change  and  7/  will  decay  in  steps  of 
approximately  cj>c !L  as  flux  quanta  enter  the  loop.  In  the 
linear  region,  for  7/«y,  the  rate  of  change  in  the  loop 
current  is  related  to  Aw  (for  constant  Oa)  by 


dt  co 


(4) 


with  the  junction  voltage  K=/iv<J>0. 

It  is  important  to  note  from  this  discussion  that,  within 
the  constraint  of  fluxoid  quantization,  it  is  possible  for 
the  two  junctions  to  have  slightly  different  average  volt¬ 
ages  while  in  a  fixed  fluxoid  state,  phase  locked  to  the 
applied  radiation.  This  result,  contained  in  Eq.  (4),  has 
previously  been  demonstrated  with  a  modified  version  of 
the  apparatus  in  which  separate  coaxial  cables  extended 
from  each  junction  to  outside  the  cryostat.12  This  per¬ 
mitted  the  creation  of  a  small  frequency  difference  by 
Doppler  shifting  of  the  radiation  going  to  one  of  the 
junctions  and  served  to  demonstrate  that  the  sensitivity 
of  the  apparatus  is  indeed  correctly  predicted  by  Eq.  (4). 

The  sensitivity  limit  of  the  apparatus  at  a  loop  emf  of 
1  x  1 0 —22  V  over  a  typical  10-h  run  implies  that  the 
change  in  the  loop  flux  <t>  due  to  any  extraneous  sources 
must  be  less  than  2x  10-4<J>o/h.  To  provide  shielding 
against  such  small  changes  in  the  ambient  flux  <I>a,  Mu¬ 
metal  shielding  was  used  to  reduce  the  ambient  field  to 
about  1  mG.  To  stabilize  the  remaining  flux  the  super¬ 
conducting  loop  was  enclosed  in  a  superconducting  NbTi 
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box  which  was  placed  inside  a  large  lead  can.  The  tem¬ 
perature  of  the  box  was  regulated  to  ±  50  /iK  to  prevent 
thermally  induced  flux  motion.  In  addition  both  junc¬ 
tions  were  fabricated  on  sapphire  substrates  to  minimize 
any  thermal  gradients  across  the  junction  and  prevent  a 
thermoelectric  current  from  being  induced  in  the  loop. 
Deliberate  temperature  shifts  of  1  mK  caused  changes  in 
O  of  less  than  of  lxl0_4$o  even  during  temperature 
transients. 

Additional  factors  which  could  cause  A/i  or  A/2  to 
change  and  generate  a  loop  current  //  are  changes  in  the 
microwave  power  or  frequency  or  in  the  junction  bias 
currents  /,-.  Since  there  are  common  sources  of  mi¬ 
crowave  and  bias  current  for  both  junctions,  the  changes 
mentioned  above,  except  for  the  gravitational  red  shift  in 
frequency,  will  affect  the  junctions  equally.  For  identical 
junctions  these  extraneous  effects  would  then  produce 
equal  changes  in  A/i  and  A/2  causing  no  change  in  //. 
Even  with  the  small  asymmetries  present  in  the  actual 
apparatus,  it  has  proven  possible— as  detailed  below — to 
reduce  the  effects  of  variations  in  the  bias  current  and 
microwave  sources  to  negligible  levels. 

Current  from  the  common  dc  bias  source  was  divided 
between  the  junctions  so  that  to  first  order  a  change  in 
the  source  current  produced  no  change  in  //.  The  drift  in 
the  current  source  of  1  part  in  104  during  a  run  corre¬ 
sponded  to  an  apparent  voltage  of  about  10  _23  V.  The 
frequency  drift  of  the  microwave  source,  which  was 
phase  locked  to  a  quartz-crystal  oscillator,  was  1  kHz/d. 
This  resulted  in  a  coherent  drift  in  the  voltages  of  the 
two  Josephson  batteries  of  about  a  picovolt  over  the 
course  of  the  measurement.  However,  the  symmetric 
configuration,  with  equal  microwave  path  lengths  to  the 
two  junctions,  reduced  the  flux  change  through  the  loop 
because  of  this  coherent  drift  to  a  rate  of  4.2  x  10_5Oo/ 
d,  corresponding  to  an  apparent  loop  emf  of  less  than 
10 _24  V.  The  power  output  of  the  microwave  source 
was  regulated  to  about  2  parts  in  103  by  use  of  a  leveling 
loop.  However,  the  power  reaching  the  junctions  varied 
by  about  1.5%  during  the  measurement  because  of  a 
change  in  the  attenuation  of  the  coaxial  cable  as  the  heli¬ 
um  level  in  the  cryostat  dropped.  This  variation,  which 
was  the  largest  source  of  sample-related  flux  drift,  pro¬ 
duced  an  estimated  change  in  O  of  5x  10  — 4<t»o  over  12  h 
to  give  an  apparent  loop  emf  of  0.4x  10  _22  V. 

The  loop  flux  was  monitored  for  continuous  10-h 
periods  with  use  of  a  commercial  rf  SQUID  magnetome¬ 
ter  located  in  the  helium  bath.  The  magnetometer  out¬ 
put  d>m  (referred  to  the  sample  loop)  versus  time  is 
shown  in  Fig.  2,  curve  a,  for  the  loop  junctions  phase 
locked  on  their  ninth-order  steps  to  16.8-GHz  radiation. 
This  gives  a  junction  voltage  of  F— 300/iV,  and  a 
difference  between  the  two  junction  voltages,  as  a  result 
of  the  gravitational  red  shift,  of  2.35xl0-21  V.  <J>m  is 
essentially  time  independent  except  near  the  end  of  the 
run  when  the  helium  level  is  near  the  rf  SQUID.  This 


FIG.  2.  Curve  a,  vs  time  with  the  junctions  phase  locked 
to  16.8-GHz_  radiation  biased  at  300  juV  on  the  ninth-order 
step. '  Curve  b,  <t>m  vs  time  with  the  junctions  biased  on  the 
zeroth  step.  Curve  c,  ^(F^SOO  pV)  —  4>m(F“0 pV)  vs 
time.  The  dashed  line  shows  the  expected  signal  if  jr  rather 
than  p  remained  constant  along  the  superconducting  wire.  The 
solid  lines  in  a  and  c  are  the  respective  best  linear  fits  as  de¬ 
scribed  in  the  text. 


dependence  of  on  time  is  highly  reproducible  from 
run  to  run.  To  check  that  the  residual  signal  is  due  to 
the  rf  SQUID,  the  measurement  was  repeated  with  the 
junction  bias  current  reduced  so  that  the  junctions  were 
on  their  zeroth  step  (i.e.,  F  — 0  /rV).  The  resulting  vari¬ 
ation  of  <t>m  with  time,  shown  in  Fig.  2,  curve  b,  is  essen¬ 
tially  identical  to  that  in  Fig.  2,  curve  a,  indicating  that 
there  is  no  measurable  change  in  the  flux  through  the 
sample  loop  during  the  run.  The  maximum  flux  change 
through  the  sample  loop  as  a  result  of  relativistic  effects 
is  just  the  difference  between  curves  a  and  b,  i.e., 
d>n,(F  =  300  /iV) -OmfF^O  /iV)  as  shown  in  Fig.  2, 
curve  c. 

A  linear-regression  analysis  of  the  data  in  Fig.  2,  curve 
c,  taking  into  account  the  presence  of  1//  noise  in  the 
SQUID  output,  gives  for  the  magnitude  and  error  limits 
(2cr)  of  the  loop  voltage  (3±6)xl0-23  V.  Nearly  the 
same  limits  are  obtained  from  the  data  in  Fig.  2,  curve  a, 
if  the  final  4  h  of  data  at  low  helium  level  are  discarded. 
In  this  case  we  obtain  (1  ±8)x  10  ~23  V. 

We  thus  conclude  that  the  net  emf  in  a  superconduct¬ 
ing  loop  containing  two  Josephson-effcct  batteries  at 
different  gravitational  potentials  is  not  the  difference  in 
the  battery  electrochemical  potentials,  AF=XF,  but 
rather  is  indistinguishable  from  zero,  consistent  with  the 
relativistic  prediction  that  fi  rather  than  /r  is  constant 
along  the  superconducting  wires.  In  particular  we  have 
tested  the  hypothesis  that  the  potential  difference  be¬ 
tween  the  wires  varies  with  height  as  F(z)“F(0) 
x  (1  —  X9).  Our  results  imply  that  X/X9  =  1  ±  0.04,  as 
compared  with  the  predicted  equality  of  X  and  X, 
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Measurements  have  been  made  of  the  transition  rates  between  fluxoid  states  of  a  system  with  two 
Josephson  junctions  having  two  macroscopic  degrees  of  freedom.  The  data,  taken  over  a  temperature 
range  of  0.85  to  4.41  K,  are  in  excellent  agreement  with  the  predictions  for  thermally  activated  transi¬ 
tions.  No  evidence  is  seen  for  recently  reported  apparent  rate  suppression  below  that  predicted  for 
thermal  activation. 


PACS  numbers:  74.50.+r,  05.40.+j,  85.2S.Cp 


Josephson  junctions  have  proven  to  be  excellent  sys¬ 
tems  in  which  to  study  the  dynamics,  both  classical  and 
quantum,  of  a  macroscopic  coordinate  in  a  metastable 
potential.  In  general,  careful  experiments  have  yielded 
excellent  agreement  with  theory  for  systems  with  one  de¬ 
gree  of  freedom,  i.e.,  a  single,  small  junction.  Recent  re¬ 
sults  by  Sharifi,  Gavilano,  and  Van  Harlingen'  7  (SGV) 
on  thermal  activation  from  the  zero-voltage  state  of  a  dc 
SQUID  (a  system  with  two  degrees  of  freedom)  have, 
hov/ever,  shown  marked  disagreement  with  classical 
theory  for  two-dimensional  (2D)  potentials  and  have 
been  interpreted  as  showing  a  dramatic  suppression  of 
the  thermal  activation  rate  from  that  predicted.  In  this 
Letter  we  report  the  results  of  thermal  activation  mea¬ 
surements  on  a  similar  two-dimensional  system  of 
Josephson  junctions  (a  combined  rf-dc  SQUID)  in  which 
the  potential  can  be  independently  determined  to  very 
high  accuracy. 

The  system  studied  here,  as  shown  in  Fig.  1,  consists 
of  two  Josephson  junctions  in  parallel  in  a  low  induc¬ 
tance  (/)  superconducting  loop  (the  dc  SQUID)  which  is 
incorporated  as  part  of  a  larger  inductance  (L)  super¬ 
conducting  loop  (the  rf  SQUID).  The  macroscopic 
dynamical  variables  which  describe  this  system  are  the 
fluxes  4>  (which  can  be  monitored)  and  Ode  linking,  re¬ 
spectively,  the  rf  and  dc  SQUID’s.  The  fluxes  applied  to 
these  loops,  0*  and  O^c,  can  be  independently  con¬ 
trolled.  The  dynamics  of  this  system  can  be  visualized  as 
motion  in  a  2D  potential,  which  for  symmetric  dc 
SQUID  having  equal  junctions  is  given  by 

U(<p,<p &) “--r-l  t  (<P-<PX)2+  2  yW-V’x dc)2 
4  n  L 

PcOSjPdc  COS<?]  ,  (1) 


3  0 


Magneto¬ 

l _ 

meter 

FIG.  1.  Schematic  of  SQUID  sample  and  coupling  coils. 


with  a  damping  determined  both  by  the  junctions'  resis¬ 
tance  and  the  environment.  Here,  Oo  is  the  flux  quan¬ 
tum,  y=L/l  is  the  ratio  of  the  loop  inductances,  and 
P=2jtLIJ<bo>  where  Ic  is  the  sum  of  the  junctions’  criti¬ 
cal  currents.  V>o“*(I)u[2ff/<X>o].  Figure  2  shows  a  plot  of 
this  potential  for  >3  “>4.14  and  y  “19.5  (the  case  for  the 
junctions  measured  here)  with  O*  “0.65O0  and  4>^dc“0. 
For  these  parameters  the  potential  consists  of  two  local 
minima  connected  by  a  saddle  point.  For  the  values  of  p 
and  y  here,  there  are  at  most  two  local  minima  with 
C/(<3>,« I>dC)  increasing  roughly  parabolically  at  large  O 
and  d>dc.  For  values  of  0*dc  >  0  the  barrier  amplitude  is 
reduced,  finally  going  to  zero  for  0*dc  ”C.4050o  (for 


FIG.  2.  The  2D  potential  U (<&,<$>&)  of  our  sample  SQUID 
with  symmetric  dc  SQUID  having  identical  junctions.  The 
mutual  inductance  between  two  loops  is  taken  to  be  1/2.  The 
potential  for  p  "4. 14,  y—  1 9.5,  O*  "O.650o,  and  0Xdc  "*0  show¬ 
ing  (a)  equipotential  contours  at  the  above  value  of  external 
fluxes  and  (b)  L/(<X>,0).  U  continues  to  increase  monotonically 
(roughly  parabolically)  for  greater  values  of  O  and  <t>dc  than 
those  shown. 
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id>o)  so  that  no  metastable  states  exist.  Also  for 
^xdc  >  0  the  path  linking  the  wells  through  the  saddle 
point  is  no  longer  a  straight  line.  The  details  of  this  po¬ 
tential  along  with  the  techniques  used  for  the  measure¬ 
ment  of  its  parameters  will  be  presented  in  a  forthcom¬ 
ing  extended  paper. 

The  longitudinal  (transverse)  small  oscillation  fre¬ 
quencies  parallel  (perpendicular)  to  the  path  from  the 
well  through  the  saddle  point  are  oiw  (colw)  in  the  well 
and  6>|j  (wts)  at  the  saddle  point.  The  predicted  rate  of 
thermally  activated  escape  from  such  a  two-dimensional 
well  is  given  in  the  classical  limit  by  Ben-Jacob  et  al. 3  to 
be 

r2D-^exp[— At/(p,¥>dc)Ai?7’],  (2) 

where  AU  is  the  energy  difference  between  the  metasta¬ 
ble  potential  well  and  the  saddle  point  and  ft***fto a,. 
ft  ‘m(oiwo)tJ(o,s  and  a,  is  a  damping-dependent  suppres¬ 
sion  of  the  rate  which  can  be  much  less  than  one  in  both 
high-  and  low-damping  limits,  i.e.,  for  a>  much  greater  or 
much  less  than  RC.  R  and  C  are,  respectively,  the 
junction’s  resistance  and  capacitance.  For  moderate  to 
heavy  damping  a,  “*(l+<72/4)  —  <7/2,  where  G  — 1/ 
RCa>is.  In  the  extreme  underdamped  limit,  a, « 1  is  also 
predicted  since 

AU 
kBT 

Recent  results  for  single-current  biased  junctions  provide 
evidence  for  this  extreme  underdamped  limit  with  a 
damping  ^ccexp(— AU/kBT),  where  AU  is  the  super¬ 
conducting  energy  gap,  as  expected  for  quasiparticles.4'6 
However,  the  extension  of  these  results  to  multidimen¬ 
sional  potentials  is  still  open  to  question.7 

All  measurements  were  made  in  highly  shielded  cryo¬ 
stats  (either  4He  or  dilution  refrigerator).  All  signals 
were  inductively  coupled  to  the  sample,  which  was  en¬ 
closed  in  a  NbTi  can,  through  leads  having  low- 
temperature  low-pass  filters.  For  low-temperature  mea¬ 
surements  a  resistive  shunt,  heat  sunk  to  the  dilution  re¬ 
frigerator  mixing  chamber,  was  placed  across  the  mag¬ 
netometer  input  to  filter  any  high-frequency  signals  gen¬ 
erated  by  the  SQUID  magnetometer.  The  junctions 
used  were  nominal  1  x  l  ^m2  Nb/AhOj/Nb  tunnel  junc¬ 
tions  with  very  low  subgap  leakage.  The  I-V  curve  of  an 
unshunted  (without  L )  SQUID  coprocessed  with  the 
sample  SQUID  is  shown  ;n  the  inset  of  Fig.  3(b).  These 
data  give  a  single-junction  differential  resistance  of 
Rn  “320  ft  above  the  gap  and  a  critical  current  of 
Iem 5.8  /iA.  From  the  dc  SQUID  resonance  of  this 
unshunted  SQUID  with  d>xdC>0,  a  capacitance  of 
C**46±6  fF  was  determined  in  agreement  with  the 
specific  capacitance  obtained  by  Lichtenberger  et  al.  for 
similar  junctions.8  The  inductance  /  of  the  dc  SQUID 
was  determined  from  the  sample  as  discussed  below. 


2nwu  i 
RUcoiw 


0  12  3  4  5 

T  (K) 


FIG.  3.  Thermal  activation  rates  between  fluxoid  states  for 
the  symmetric  potential  (<t>x  •“  ?  <t>o)  shown  as  the  inset  in  (a), 
(a)  Rate  vs  l/T  for  fixed  barrier  height  with  best  linear  fit 
(line),  (b)  Barrier  height  vs  T  for  fixed  rate  (T-l  s-1).  For 
points  indicated  by  squares,  data  have  been  corrected  for  tem¬ 
perature  dependence  of  critical  current.  Best  fit  by  Eq.  (2)  as¬ 
suming  damping  due  to  R„  (solid  line),  or  damping  proportion¬ 
al  to  exp[-A U/keT]  from  quasiparticles  (dashed  line).  The 
inset  of  (b)  shows  the  dc  SQUID  I-V  curve  for  Tm  1.5  K, 
d>xdc™0  indicating  low  subgap  leakage. 


For  the  special  case  where  <t>x  is  a  half-integer  flux 
quantum,  the  potential  is  symmetric  about  so 

that  At/(fl>xdc)  is  the  same  for  each  well  [see  inset,  Fig. 
3(a)],  varying  from  a  maximum  of  143  K  for  d>xdc*-0  to 
0  for  <J>xdc  >  0.405<t>o,  where  the  potential  has  only  one 
minimum.  Figure  3(a)  shows  the  measured  rate  at 
which  the  sample  hops  between  two  equal  energy  fluxoid 
states  versus  temperature  with  constant  AU  (fixed  d>xdc). 
If  ft  is  temperature  independent  in  this  range,  one  ex¬ 
pects  from  Eq.  (2)  that  a  linear  relationship  exists  be¬ 
tween  ln(r)  and  l/T  since 

ln(r)  ^AU/kT  -Hn(ft/2/r) .  (3) 

As  Fig.  3  shows,  the  observed  dependence,  is  in  fact 
linear,  giving  best-fit  values  of  AU  •*•32.2  ±0.4  K  and 
ln(ft/2?r)  "*23.9  ±0.3.  These  agree  well  with  the  values 
of  AU  “  33.5  ±  3  K  and  In  ( ft/ 2n)  24.4  ±0.2  calculated 
from  measured  sample  parameters,  as  discussed  below. 
Here  almml  has  been  used;  taking  the  damping  to  be 
I /Rn  would  give  a,  =  0.9— an  insignificant  correction. 

The  greatest  uncertainty  in  the  calculated  value  of  AU 
comes  from  the  knowledge  of  L,  which  has  been  calculat¬ 
ed  from  the  loop  geometry,  giving  Z.“"230  pH  ±10%. 
Measurements  on  single-junction  SQUID’s9,10  have 
shown  such  calculations  to  be  reliable  for  the  geometry 
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used  here.  Both  y  and  the  difference  in  the  junctions’ 
critical  currents,  S=(J2~ /|)//c,  can  be  determined  from 
the  equilibrium  dependence  of  <t>  on  O*  measured  when 
d>,dc  is  a  half-integer  flux  quantum  so  that  the  system 
has  no  metastable  states.  These  data  give  y—19.5  and 
S  <  5%  justifying  the  treatment  of  the  junction  parame¬ 
ters  as  equal,  p  is  determined  by  measuring  the  value  of 
<t>x  at  which  the  system  escapes  from  its  metastable  state 
for  d>xdc“0.  Since  this  escape  will  occur  slightly  before 
the  barrier  is  reduced  to  zero,  a  small,  model-dependent 
correction  to  p  based  on  Eq.  (2)  is  required.  This  correc¬ 
tion,  at  most,  increases  the  calculated  AU  above  by  5% 
over  that  which  would  be  obtained  by  assuming  the  es¬ 
cape  occurred  for  AU  — 0. 

The  transition  rate  as  shown  in  Fig.  3(a)  is  only 
measurable  over  a  narrow  range  of  temperature  for  fixed 
AU.  To  measure  the  temperature  dependence  of  this  rate 
over  a  broad  range  of  temperature,  AU  was  adjusted  at 
each  temperature  to  give  approximately  the  same  rate, 
10  <  T  <  50.  These  data  are  plotted  as  AU'  vs  T  in  Fig. 
3(b).  Here  AJ/'==A{/+fcfi7Tn[r/l  s-1]  is  the  measured 
barrier  normalized  for  the  small  variations  in  ln(r) 
among  the  points.  Again  <PX  -  |  <J>0  and  A£/(<I>,x(jc)  was 
adjusted  by  varying  and  calculated  from  the  param¬ 
eters  as  above. 

The  theoretical  curve  for  thermal  activation  has  been 
plotted  using  Eq.  (3)  taking  the  damping  to  be  1  /R„ 
which  gives  a,  ^  0.9.  The  uncertainty  in  the  indepen¬ 
dently  measured  values  of  AU  ±  10%  implies  an  accept¬ 
able  fit  to  the  data  in  Fig.  3(b)  for  any  temperature- 
independent  damping  less  than  0.05  ft  _l  if  the  moderate 
damping  expression  for  a,  is  used.  The  dashed  line 
in  Fig.  3(b)  is  calculated  assuming  the  damping  17 
«exp(— A U/ksT)  as  expected  for  quasiparticles  and 
using  the  extreme  underdamped  expression  for  a,.  This 
is  not  consistent  with  our  data,  however,  we  emphasize 
that  this  result  may  imply  only  that  extrinsic  damping 
(e.g.,  losses  from  the  SQUID  loop)  dominates. 

We  conclude  that  the  data  presented  in  Fig.  3  arc  en¬ 
tirely  consistent  with  thermal  activation  [Eq.  (2)]  and,  in 
particular,  do  not  admit  an  effective  barrier  height 

TABLE  1.  Characteristic  energies  (in  units  of  Kelvin)  of  the 
symmetric  double-well  potential  ($>.,  —  [  <t\)  which  describes 
the  sjstcm  for  data  in  Fig.  3.  At  each  temperature  has 
been  adjusted  to  give  the  values  of  AU  shown. 


7- 0.85  K 

7-4.41  K 

AU 

18 

92 

han. 

1.85 

3.15 

h<Utw 

14.6 

14.8 

hcou 

1.36 

2.61 

ha>u 

14.4 

14.2 

AUlhmu 

9.7 

29 

AU/ho},m(ii) 

-1.2 

-6 

FIG.  4.  Mean  (<&*)  (circles)  and  width  a  (squares)  of  the 
thermal  activation  probability  distribution  (see  text)  vs  T  at 
Ocdc”  0.1224>o.  The  solid  and  dashed  lines  are,  respectively, 
the  predicted  mean  and  width  from  Eq.  (2)  using  parameters 
given  in  text.  Inset.  A  typical  distribution  used  to  calculate 
(4>x)  and  a  at  one  temperature. 

AU*  — 2.5AU  as  needed  by  SGV  to  explain  their  data  in 
the  thermal  activation  regime.  The  energies  associated 
with  the  symmetric  potential  over  the  temperature  range 
of  the  measurements  are  summarized  in  Table  I.  For 
comparison,  we  note  that  both  for  these  data  and  for 
those  of  SGV,  hcoi^kaT,  whereas  hm^ksT  with 
<o,la>i «  8.  Even  for  our  lowest-temperature  data  there 
are  at  least  ten  energy  levels  in  the  wells,  satisfying  the 
validity  criteria  for  the  classical  [Eq.  (2)]  and  semiclassi- 
cal"  rate  calculations.  The  expected  crossover  tempera¬ 
ture  Tx  for  tunneling-dominated  transitions  is  Tx  -0.35 
K  in  the  zero-damping  limit.  Thus  quantum  correction11 
to  Eq.  (2)  should  be  negligible  for  our  data.  Naor, 
Tesche,  and  Ketchen12  have  also  reported  excellent 
agreement  witi.  thermal  activation  theory  for  their  mea¬ 
surements  for  transitions  for  r>  1.7  K  in  a  highly 
damped,  current-biased  dc  SQUID  with  <ot/<oi  ~  1. 

A  more  direct  comparison  with  the  results  of  SGV  is 
provided  by  the  data  in  Fig.  4.  Here,  for  Oj„ic— 0.1220o, 
Ox  is  increased  at  a  constant  rate  reducing  the  barrier 
height  until  the  SQUID  undergoes  a  change  of  fluxoid 
state  by  escaping  from  a  metastable  well.  Measurements 
of  this  type  on  a  single-junction  SQUID  have  been  de¬ 
scribed  in  detail  elsewhere.9 10  These  measurements  are 
repeated  approximately  104  times  at  a  fixed  temperature 
giving  a  histogram  of  transition  probability  versus  O* 
(see  Fig.  4,  inset).  The  first  and  second  moments  of  this 
distribution,  (<t>*>  and  <7,  respectively,  are  predicted  to 
vary  with  temperature13  as  shown  by  the  solid  and 
dashed  lines  in  Fig.  4  assuming  a  potential  and  transition 
rate  given  by  Eqs.  (1)  and  (2)  and  using  the  sample  pa¬ 
rameters  above. 

In  calculating  (O*)  vs  T  the  critical  applied  flux  Oxc  at 
which  A£/— 0  is  used  as  a  temperature-independent 
fitting  parameter,  Ow  -0.8820o,  which  is  consistent 
with  P  measured  for  Oxdc— 0.  As  can  be  seen,  the  agree- 
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ment  with  thermal  activation  predictions  is  excellent. 
This  contrasts  with  similar  measurements  by  SGV  in 
which  they  found  for  their  system  that  d(<t>x)/dT  was 
about  a  factor  of  3  less  than  predicted.  The  width  cr  also 
varies  with  temperature  as  expected,  although  it  is  about 
5%  greater  than  predicted.  This  is  consistent  with  the 
excess  width  expected  from  the  Nyquit  flux  noise  on  the 
sample  due  to  the  shunt  resistor  across  the  magnetome¬ 
ter.  The  temperature  of  this  shunt  varied  from  50  to  180 
mK  as  the  sample  temperature  was  increased. 

Measurements  have  been  made  of  thermally  activated 
transitions  between  fluxoid  states  of  a  SQUID  system 
having  two  macroscopic  degrees  of  freedom,  in  which  the 
energy  potentials  could  be  accurately  determined  and 
where  existing  classical  and  semiclassical  theory  should 
be  valid.  Excellent  agreement  has  been  found  with  the 
prediction  of  thermal  activation  theory.  In  particular,  no 
hint  was  seen  of  the  very  large  apparent  rate  suppression 
recently  reported  in  a  similar  system  by  SGV.  In  view  of 
the  similarities  of  the  potentials  for  these  two  measure¬ 
ments  (in  particular,  the  ratios  of  the  longitudinal  and 
transverse  level  spacing  to  kBT),  our  results  should  serve 
to  focus  on  the  key  features  of  the  system  used  by  SGV 
responsible  for  their  anomalous  results.  lor  example, 
our  results  indicate  that  such  anomalies  cannot  be  a  gen¬ 
eral  feature  of  two-dimensional  potentials  even  in  the 
limit  where  ha>,^>kBT  so  that  level  quantization  of  the 
transverse  mode  is  important. 
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